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ABSTRACT 


Assuming that cosmic ray deuterons and helium-3 nuclei are 
of secondary origin, we show that a unique determination of both 
the cosmic ray path length and the residual interplanetary field 
modulation at solar minimum may be made from a comparison of the 
calculated and measured intensities of these two nuclei. This 
determination does not depend on any assumptions regarding either 
the source spectra or the unmodulated proton-to-alpha particle 
ratio of the primary cosmic rays. The production of deuterium 
and helium-3 by cosmic ray interactions in the galaxy is calcula- 
ted, considering energy-dependent cross-sections, interaction 
kinematics and demodulated cosmic ray spectra. The resulting 
flux at the earth is obtained by taking into account leakage 
from the galaxy, ionization losses, nuclear breakup and modu- 
lation. From a comparison of these calculations with the meas- 
ured deuterium and helium-3 intensities at the earth, we conclude 
that within the experimehtal uncertainties all the data can be 
understood in terms of an energy-independent cosmic ray path- 
length of 4 + 1 g/cm^ and a residual interplanetary field modu- 
lation at solar minimum which above 600 MVlSofthe form exp(-T]/Rp) 
with r) = 0.35 + 0.15 BV , where R and g are rigidity and velocity 
respectively . 


COSMIC RAY DEUTERIUM AND HELIUM-3 OF SECONDARY ORIGIN 
AND THE RESIDUAL MODULATION OF COSMIC RAYS 

I. INTRODUCTION 

The abundance of deuterium and helium-3 in the cosmic 
rays was found to be significantly larger than that expected 
from the study of the universal abundance of elements. It 
has been suggested, therefore, that these isotopes are of 
secondary origin, produced by nuclear interactions of cosmic 
rays with the interstellar gas and that the comparison of the 
observed and calculated intensities of deuterium and helium-3 
could determine the mean amount of material traversed by cosmic 
rays in interstellar space. 

O 

Previous studies of d and He production were made by 
Hayakawa et al (1958) , Foster and Mulvey (1963) , Badhwar and 
Daniel (1963), Dahanayake et al (1964), Kuzhevskii (1966), 
Fichtel and Reames (1966), Badhwar and Kaplon (1966), Biswas 
et al (1966) , Meyer et al (1968) , and Biswas et al (1968) . 

The mean cosmic ray path lengths deduced from these studies 
range from a few g/cm to more than 10 g/cm^ and they may or 
may not be energy dependent. Even though these values are of 
the same order as the path length deduced from the observed 
abundances of Li, Be and B (Shapiro and Silbefll*^ 1968) , the 


- 2 - 


determinat ion of the amount of material traversed by cosmic 
rays at low energies is seriously complicated by an unknown 
residual solar modulation. 

The existance of such a residual modulation, active even 
at solar minimum, is indicated by the measurements of an inter- 
planetary gradient, between the orbits of the earth and Mars, 
in the proton and alpha particle fluxes (O'Gallagher and Simpson, 
1967; 0’ Gallagher, 1967). The measurement of this gradient, 
however, determines only the relative modulation over a distance 
of roughly an astronomical unit, and not the total residual mod- 
ulation, which short of direct measurements over large distances 
in the solar system, remains essentially undetermined. 

In a preliminary report of rather extensive calculations 
of deuterium and helium-3 production (Ramaty and Lingenfelter , 
1968), we suggested that because of the secondary origin of these 
nuclei and because of their different charge-to-mass ratio, a 
comparative study of d and He^ in the cosmic rays allows an es- 
sentially independent determination of both the amount of ma- 
terial traversed by cosmic rays and their total residual modu- 
lation in the interplanetary magnetic field. Unlike previous 
attempts to estimate the total modulation, this technique re- 
lies solely on the assumption that the observed d and He^ are 
of secondary origin, and does not depend on any arbitrary assump- 
tion about the charge ratios or spectral shapes of the primary 
cosmic rays either in interstellar space or at their sources. 
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In the present paper we give a detailed discussion of these 
calculations and conclusions. 

The production spectra of deuterium and helium-3, re- 
sulting from the interaction of cosmic ray protons, alpha — 
particles and medium nuclei with the interstellar gas, are 
calculated by taking into account the energy dependence of 
the cross-sections, the kinematics and the secondary energy 
distributions of the various d and He^ producing interactions. 
Using these spectra, the equilibrium d and He^ fluxes are eval- 
uated by considering leakage from the galaxy, ionization losses 
and nuclear breakup. The residual solar modulation is taken 
into account by first demodulating the observed primary cosmic 
ray fluxes at the earth, thenky using these demodulated inten- 
sities to calculate the production of d and He^ in interstellar 
space and finally by modulating the resultant secondary equili- 
brium fluxes with the same modulating function as used in the 
initial demodulation process. By comparing the resultant d 
and He^ spectra with the observations, it is demonstrated that 
the measured deuteron and helium-3 fluxes at the earth can be 
best understood in terms of an energy-independent path length 
and a residual solar modulation which is both velocity and ri- 
gidity dependent. 



II. DEUTERON AND HELIUM-3 PRODUCTION 


The rate of production of secondary nuclei by cosmic ray 
interactions in the interstellar material may be written as 

1 (Q = JjEi l (E)*U«:WF l (E,E f ) (1) 

l 

where q g is the production spectrum of secondary nuclei per 
second per gram of interstellar material; E and E s are the 
energies per nucleon of the primary and secondary particles 
respectively; j is the equilibrium intensity in interstellar 
space of the interacting cosmic ray nuclei; n is the number of 
target nuclei per gram of interstellar material; <5* is the inter- 
action cross-section; and F(E,E s )dE s is the probability that a 
primary cosmic ray nucleus of energy per nucleon E will produce 
a secondary nucleus of energy per nucleon in dE s around E s . 

The production is summed over all interactions, i, which 
produce the secondary nucleus, s. The interactions which we 
have considered for the production of deuterium, helium-3 and 
tritium, which decays into helium-3, are listed in Table 1. 

Each of these interactions may also include the generation of 
pions . 

The interstellar material is assumed to be composed of 
hydrogen, helium and CNO nuclei in the ratio 1:10“^: 10“^ (Suess 

and Urey, 1956; Cameron, 1959) 

We shall now provide detailed discussions of the primary 
cosmic ray intensities, production cross-sections, and the 



-5- 


kinematics and secondary energy spectra to be used in the 
evaluation of equation (1) . 

Primary Cosmic Ray Intensities 

The energy spectra of the various charge components of the 
primary cosmic radiation were measured near the earth at solar 
minimum by Freier and Waddington (1965) , Balasubrahmanyan et al 
(1966 a,b), Comstock et al (1966), Fan et al (1966b), Hofmann 
and Winckler (1966), Ormes and Webber (1966) and Waddington and 
Freier (1966). The interstellar primary cosmic ray intensities 
used to evaluate equation (1) are obtained by demodulating 
these measured spectra with a modulating function, M , of the 
form {(?))> based on Parker's (1958) solar wind model. 

The modulating parameter , which is space and time dependent, 
defines the total residual solar modulation, p> is the particle 
velocity in units of c , and f(R) is a function of magnetic 
rigidity and depends on the distribution of the interplanetary 
magnetic field irregularities. It was first suggested by Dorman 
(1963) that f (R) could be approximated by f(R)=J^,for R<R 0 and 
f(R)=R for R>R #> where R c is a characteristic transition rigidity 
depending on the distribution of the magnetic irregularities. 
Subsequently, Jokipii(1966 , 1968) and Jokipii and Coleman (1968) 
demonstrated that such a form is consistent with the measured 
power spectrum of the interplanetary field variations, with 
R Q essent ially time dependent and of the order of several 
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hundred MV in 1965. With this functional form for f (R) , the 
modulating function mentioned above can be written as 

M = exp (-•£ / , R<J?« 

( 2 ) 

ex F (“ \ ! R > 

Gloeckler and Jokipii (1966) , Balasubrahmanyan et al (1967) , 

O' Gallagher and O' Gallagher (1967) have also shown 

that the measured temporal and spatial variations of cosmic rays 
in the interplanetary medium can be explained by such a modula- 
tion, but have deduced different values for the modulating para- 
meter R Q . Since there is no general agreement as to the value 
of R q , and since the modulating parameters which define the total 
modulation may be different from the values obtained from the 
temporal and spatial variations mentioned above, we have treated 
both and R Q as free parameters to be determined from the pres- 
ent study. 

Production Cross-Sections 

The deuteron production cross-section for proton-proton 
interactions has been measured in considerable detail from 
threshold to a few Bev by Schulz (1952), Crawford (1953), 

Stevenson (1953), Durbin et al. (1951), Fields et al. (1954), 
Stadler (1954) , Baldoni et al (1962) , Guzhavin et al (1964) , 
Neganov and Parfenov (1958) , Chapman et al (1964) , Fickinger 
et al (1962), Pickup et al (1962), Smith et al (1961) and Hart 
et al (1962), and the theory of the interaction has been studied 
by Rosenfeld (1954). This cross-section data together with a 
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smoothed curve used in the present calculation are shown in 
Figure 1. 

The cross-sections for the production of deuterons, tritons 
and helium-3 nuclei from the breakup of He 4 were measured in 
proton interactions with alpha particles at 28 Mev by Wickersham 

(1957) , at 31 Mev by Bunch et al (1964) , at 32 Mev by Benveniste 

and Cork (1953) , at 53 Mev by Cairns et al (1964) , at 55 Mev by 

Hayakawa et al (1964), and at 95 Mev by Selove and Teem (1958), 

and in neutron interactions with alpha particles at 90 Mev by 
Tannenwald (1953) , and at 300 Mev by Moulthrop (1955) and Innes 
(1957). These data are shown in Figure 2. We have assumed 
that the He 4 breakup cross-sections at these energies are the 
same for incident protons as for neutrons. The measurement of 
the pickup cross-section by Tannenwald (1953) of 12+2.5 mb for 
90 Mev neutrons in the reaction rio(-»dt and by Selove and Teem 
(1958) of 15+3 mb for 95 Mev in the reaction po(-*dHe^ suggest 
that the coulomb barrier for protons is not important at these 
energies and that these cross-sections are equal. At lower 
energies, particularly near the threshold, the coulomb barrier 
is important and the proton pickup cross-section is much lower 
than that for neutrons, which shows a resonance at about 22.15 
Mev . 

Measurements of alpha particle breakup were also made at 
630 Mev by Kozadaev et al. (1960) and at 970 Mev by Riddiford 
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and Williams (1960) . These measurements differentiate between 
reactions having xwo , three or more charge particles in the 
final state but do not completely distinguish between the various 
modes. Therefore, some assumptions must be made in order 
to obtain the deuteron, triton and helium-3 yields. It was 
shown by Moulthrop (1955) and Innes (1957) that in 3-prong re- 
actions (three charged particles in the final state) the ratios 
between the alpha particle breakup modes into singly charged 
fragments, pt : pnd : dd :2p2n# 0 . 58 :0 . 26 : 0 . 12 : 0 . 04 , are the same 
in both pion and non-pion producing reactions and that these 
rates are roughly independent of energy. 

Using this relationship, and the fact that Kozadaev et al 
(1960) did measure the non-pionic pt breakup, we can separate 
their 3-prong cross-section into partial cross-sections for 
these breakup modes for both pion and non-pion producing re- 
actions. A similar treatment is applied to the data of Riddiford 
and Williams (1960) with the additional assumption that 1f°and 'W' 4 
production are roughly equal. 

In order to reduce the 2-prong inelastic reactions which 
3 3 

include the dHe and pnHe final states, we make three assump- 
tions: (a) pion production without alpha particle breakup con- 
stitutes roughly 10% of the total pion production cross-section 
at energies greater than 300 Mev , based on the measurements of 
Innes (1957); (b) Tf + and T° product ion are roughly equal; and (c) 
the dHe 3 and pnHe 3 final states are equally probable in pion 



-9- 


producing reactions, based on a similar equality found for 

♦ 

non-pionic modes in the lower energy measurements. The cross- 
t sections thus deduced are shown in parentheses in Figure 2, to- 

gether with the inferred energy-dependent partial cross-section 
for po< interactions . 

In the absence of data on helium-4 breakup in interac- 
tions, we have assumed that the same breakup modes occur and 
that their cross-sections are simply four times those in 
interactions at the same available kinetic energy in the center 
of mass frame. Since our calculations suggest that oW interac- 

Q 

tions make a non-negligible contribution to d and He production, 
a measurement of the breakup cross-sections in the energy range 
around 100 Mev would be particularly valuable. 

The deuteron production cross-section in proton interactions 
with CNO nuclei was measured at 18.5 Mev by Nadi and Riad (1964), 
at 20 Mev by Legg (1963), at 40 Mev by Kavaloski et al (1963) 
and at 190 Mev by Bailey (1956). At higher energies the cross- 
section can be estimated from the d/t yields calculated by 
Fraenkel (personal communication) and measured by Schwarzschild 
and Zupancic (1963). The tritium production cross-section for 
these interactions was measured at 43 Mev by Cherny and Pehl 
(personal communication), at 150 Mev by Brun et al (1962), at 
190 Mev by Bailey (1956). at 224, 300, 400 and 750 Mev by Honda 

> 

and Lai (1960), at 450 Mev and 2.05 Bev by Currie et al (1956), 
2.2 Bev by Fireman and Rowland (1955), and at 6.2 Bev by Currie 
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(1959) . The helium-3 production cross-section in proton-CNO 
interactions is not well measured, but the measurements of 
Bailey (1956) at 190 Mev and the Monte Carlo calculations of 

o O 

Fraenkel (personal communication) indicate that the He° /H° 
yields are of the order of 1.7 from protons of energies be- 
tween a hundred and several hundred Mev , and we have assumed 
that this ratio holds for all energies. These cross-sections, 
averaged to a mean value for CNO in the cosmic-ray nuclear 
ratio of 0.45:0.10:0.45 (Comstock et al, 1966) are shown in 
Figure 3. 

K inematics and Secondary Energy Spectra 

The distribution functions F(E,E S ) can be determined 
from measurements of the angular and energy distributions of 
deuterons and helium-3 nuclei produced in nuclear interactions. 
Consider the interaction of an incident particle of mass m^ 
and Lorentz factor, ^ = 1 + E/m^, with a stationary target 

nucleus of mass m t . Such a reaction may produce two or more 
secondary particles. Consider one of these particles and let 
p(y, , ( CoS be the probability that in the center 

of mass frame it is emitted i.ato dcos 9 * around cos Q * with 
Lorentz factor in d around ^*. * n terms of this pro- 
bability, the distribution function F ( ^ ^ ^ may be written 


as 


p(>, K * 


(3) 
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The laboratory Lorentz factor ^ is determined from the 
center of mass Lorentz factor and emission angle by 


X = x X + Xx 7 ***' 


(4) 


where 


V 


is the Lorentz factor of the center of mass 




_ +m t E ,Z + w^-vn.j 


E' 


Zm t £' 


(5) 


and E is the total available energy in the center of mass 
frame 




E' 

Substituting Jc0$©*/d| from equation (4) into equation (3) , 


we get 


= 


(IF 1 


i,* fsgl 


(7) 


If the final state consists of more than two particles, 
the Lorentz factor, 'y , of one of the secondary particles 
having mass m_ may have any value up to a maximum which is de- 
termined by the conservation of momentum and energy and is given 

by /> 

_ j t 

K n 


2m s E 


( 8 ) 


where m r is the sum of the masses of all the other secondary 
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particles. However, if the final state consists of only two 
particles , >r = ’ in '* the probability distribution 
reduces to a delta function in and equation (7) becomes 

ct v _ p( y> 

u>y - (9> 

For the purpose of evaluating the integrals over y and y* 
given in equations (1) and (3) , we have to consider the ranges 
of the incident Lorentz factor y and the secondary Lorentz fac- 
tor y* which can contribute to a given laboratory Lorentz fac- 
tor y . Since y is uniquely determined by y and vice versa, 
the permissible range of y also determines that of y . These 
ranges are determined by the requirements that | 9* I 1=. 1 , 

and that * Trt for more than two particles, or 

for only two particles in the final state respectively. 

Using e^utio*i(4), the restrictions on toS 8 imply that 

-1 S - 3-k ~ 1 


( m-' 


( 10 ) 


For a fixed y ^ this inequality restricts y and y to 
the inside of a hyperbola given by 




on 


For illustration, such hyperbolae are shown by solid lines in 
Figure 4 for 1.01, 1.1 and 2.0, corresponding to deuterons 

of roughly 10, 100 and 1000 Mev/nucleon respectively. The three 
branches of each hyperbola are separated by the points (V V» £*') 
and (> c *'i >4= Is) ' The value of CtiSQ* equals -1 on the 
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two branches which extend out to large values of y and "V y 
1* 

and C©s9 = +1 on the remaining branch. Inside of each hyper- 
bola CoS 8 has values between +1 and -1. 

The dashed curves in figure (4) are plots of "^as a function 


of 


Vc 


c » 


determined from equations (5) and (8) for the reac- 


tions H 1 (p,d)'Tr\ He 4 (p , d) He 3 and p(He 4 ,d)He 3 . At threshold 


Or* = 1 and \ = It + 

0 y « 4 c. Oct 


which is the Lorentz factor of the cen- 


ter of mass at threshold. Since at threshold E = m g + m r , 
from equation (5) we obtain 

(\nn s + >no r ) Z + 




2. rrn ± ( rvn s +^n r ) 


(12) 


As the incident energy increases above threshold, y 
and y also increase above their threshold values and "y 

0 YT> 

may or may not become greater than . The condition for 

* 

this may be seen from the asymptotic behavior of \ at large 

0 In 

incident energies, > * ~ ^ ^ • Thus if m t >m s , for ^ > 

j| & 

'y , where y ^ is the value of at which y* x y , 

The value of > can be directly evaluated from equations 
°cl 

(5) and (8) and is given by 


V - - no*) ~~ ( - rn * ) 

o c ? ~ f~7 x r ~ T. ^ TT" 


lc2 z K m s - w t) C i w r) - w s ( w? ) J (13) 

For values of ^ greater than y , secondary particles may 
be emitted in the laboratory frame into both the forward and 


backward cones. When 


t 


equals 


, particles emitted with 


& » -Y in the center of mass frame, have zero kinetic energy 
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in the laboratory frame. However, if m t <m s , y* can never 
equal y and in the laboratory frame, the secondary particle 
will always be emitted into the forward cone with nonvanishing 
energy . 

The lower limit, ^ , of the range of values of y which 
may contribute to a given ^ is determined by the first inter- 
section of y* with the hyperbola corresponding to the given ^ * 
as can be seen in Figure 4. Depending on whether y > or 
X < ’ this corresponds to particles being emitted in the 

center of mass frame into the forward or backward directions. 

The upper limit on y , however, is not necessarily finite. 
The asymptotic forms of the two branches of the hyperbola cor- 
responding to a given y are and, as men- 
tioned above, for large y , y* ^ ^ There- 

fore, if 


x _ [y 2 -- 1 < an* < x + ry r TT 

•c v*s ic V 


(14) 


will remain always inside the hyperbola corresponding to 
y . This inequality is equivalent to y$> where 

is equal to £(m t /m s + m s /m t ) and at large incident energies 
is the asymptotic value of the Lorentz factor of secondary par- 
ticles, y , emitted in the backward direction in the center of 
mass frame. Because of the relativistic addition of velocities, 
y Qt) is independent of the incident energy and therefore pro- 
jectiles of the greatest energies may always produce secondary 



15 


particles with yy* ), provided they are emitted close 

enough to the backward direction in the center of mass frame. 

In these cases, the upper limit, y cu > becomes infinite. 

When \ S K ^ S ( 00 ) we have to consider separately the cases 

in which m > m or in < m . For in > m and two particles 
t s t^s t s 

in the final state, the upper limit is determined by the 

second intersection of y * and the hyperbola corresponding 

q m 

to . If there are more than two particles in the final 


state, however, 


becomes infinite. On the other hand, for 


in £ m , -v, is determined by the second intersection of \ * 
t ^ s Jcu 0 m 

with the hyperbola corresponding to regardless of the num- 

ber of particles in the final state. 

Finally, a further distinction exists between reactions 

having m. greater or smaller than m , namely that for m > m , 
t S X s 

the curve 'V * will intersect every hyperbola at least once, 

0m 

whereas for m < m , y has to be greater than a minimal value, 
X S AS 

^ s (min), in order that the corresponding hyperbola be inter- 
sected by -V* . This can be seen by composing the curves of 

4 3 4 3 p* 

for the reactions H € (p, d)He" and p(He ,d)He in figure 
This minimal laboratory Lorentz factor is obtained when 
and secondary particle is emitted into the 

backward direction in the center of mass frame , Q * = -IT . 
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From equation (4) with ^, g * - 


* and 0* = -'B', we find 


that the derivative of 


h 


vanishes when 


J'm V nt ~ » 


l fc* - 1 ^ j > c i^ r ' 


- o 


( 15 ) 


The first term will never vanish since '^ m *< ^ c > the laboratory 

Lorentz factor \ will go through a minimum when the second term 

0 s 

vanishes. The values of \ and \ * which satisfy equation 

f c 

(15), (min) and (min), are obtained from equations (5) 

and (8) and are given by 


A / \ Wi ( + ~ ~ ~ l 

f c l H "V ~ ~ r~ ■ ■ — ~ r 

2 W) -t\j ^ W s z - W,?) ~ 


(16) 


and 


A *, w r ( w> t 2 -+vy> r 2 '-.»o/'- w>) h, r *) 

V (m,V) * „.=■■ - ===== — = = 

** lrm S / - rfi T~) - ( /ry, b 3 ~- ™o~) 


(17) 


The incident kinetic energies per nucleon E 


th 


E 


mm 


corresponding to ^ cz and y c ( min )> aS wel1 as the secondar y 

laboratory kinetic energies per nucleon E (th) , E (min) and 

s s 

E g («>), corresponding to y ct > y s (min) and ^(00) are shown in 


Table 2. (Note that at threshold 


J's fct 


). 
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< 


For illustration, the range of possible secondary 

energies as a function of incident energy is shown in Figure 

( 5 ) for a variety of deuteron-producing reactions. The curves 

represent the kinematic limits determined from equation (4) 

* 

for cos 0 = i 1 and form the envelope of permissible secondary 

* 

energies corresponding to intermediate values of cos© . The 

significance of the various minima and asymptotic limits listed 

in Table 2 may also be seen in Figure (5) . 

The center of mass angular distribution of deuterons from 

the reaction H(p, cl )1T^ was well measured and has been studied 

3 3 

by Rosenfeld (1954) who showed that the (— , 7 ^) interaction was 

dominant and that the distribution was therefore given by 
1 2 * 

- + cos 0 . Since this reaction also results in two bodies 

in the final state, the deuteron distribution function is then 
given directly from the angular distribution and equation (9). 


- 


/ z * 

\/5 4 Cos Q 


(li) 




\ 


The center of mass angular distribution of deuterons from 
4 3 

the reaction He (p,d)He was measured at 31 Mev by Bunch et 
al . (1964), at 55 Mev by Hayakawa et al. (1964), and at 95 Mev 

by Selove and Teem (1958). These measurements are shown in 
Figure (6) and they can be fitted in detail by the model of Smith 
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and Ivash (1962) based on the distorted wave Born approxi- 
mation with diffuse-well nuclear optical potentials. For 
the purposes of this calculation, however, we have fitted 
the distributions by functions of the form: 


P(\C0id') = e *pE ' 0 ~ — - 9 V L 

C °S0, f | - exp(- 2 /cos &o)] 


(19) 


where cosQ^ depends on the incident Lorentz factor^. In 
the incident energy region of 30 to 90 Mev/nucleon, cos Q 0 
can be well fitted by 

Cos 0 O - exp[ (2°) 


and we assumed this energy dependence for all values of 
considered . 


t 


Since this reaction results in two bodies in the final 

state, the distribution function F( Ws ) for both deuterons 

and helium-3 nuclei is determined from equations (9) and (19) . 

4 3 

The parameter £ is +1 in the reactions He (p,d)He and 

1 3 I 13 

H ( o( , He )d and is -1 in the reactions H ( o( , d)He and 

He 4 (p,He 3 )d. 

Above the pion production threshold the more complex 

3 o -f 

reactions leading to dHe ir and dtir also become signifi- 
cant', as can be seen from the cross-sections shown in figure ( 2 ). 
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It has been shown by Riddiford and Williams (1960), however, 

that pions produced in proton-helium reactions are ejected 

as from the collision of free nucleons. Because of this, we 

3 o -f“ 

have assumed that, kinematically, the dHe IT and dtir re- 
actions are similar to the reaction H(p, <J >*, and therefore 
the distribution functions F(^,,^ s ) are given by equation (18), 
where and are determined by equations (5) and (8) with 

m = m = m,m = m , and m = m 
i t p s d r 'ir 

The energy spectra in the laboratory frame of deuterons 

4 4 

from pnd breakup of He and of tritons from pt breakup of He 

were measured by Tannenwald (1953) and by Innes (1957) for 

incident neutrons of 90 and 300 Mev, respectively. These data 

are shown in Figure (7) j normalized to a unit integral. As can 

be seen, within the experimental uncertainties, the triton 

4 

spectrum from the multibody pt breakup of He is independent 
of the incident neutron energy and can be fitted by a simple 
exponential of the form: exp^-E^/E where is the triton 

kinetic energy per nucleon and 1 /E q is equal to 0.27 (Mev/ 
nucleon) - '*'. Similarly, the deuteron spectrum may be fitted 
bj** exponential with 1/E = : 0.063 (Mev/nucleon) ■*". 

These spectra indicate that in the laboratory frame the 
multinucleon fragments tend to be produced predominantly with 
low energies. This will happen if the angular and Lorentz- 
factor distributions in the center of mass are peaked in 
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the backward direction and around Lorentz factors which are 
about equal to the Lorentz factor of the center of mass. As- 
suming that the characteristic energy, E 0 , is independent of 
the incident energy, the exponential distribution in kinetic 
energy per nucleon may be transformed into an exponential dis- 
tribution in the center of mass angle, 0* . We assumed, there- 
fore. that for the multibody reactions, He 4 (p , p*) npd ,dd , pt and 
nHe 3 , the probability distribution function y* > COS G*') 

is given by 


) Cos e„ [ | - exp ( -2/cos6’ 0 '> ] Us *S / ( 21 ) 

where in po( interactions, the parameter £ equals -1 and ^ s ^ 


( 21 ) 


Wc 


Y s * for ) and in °<p 


interactions, £ = +1, and the mean Lorentz factor A* 

is the same as in po< interactions for the same total avail- 


able energy in the center of mass. 


if 

Using equations (4) and (9) with both and re- 
placed by , the laboratory distribution function 


becomes 




\ COS©, 


V 


( 22 ) 


' - exj> ( - 2 / COS ffp-i ' ] 

In terms of the characteristic energy ; E Q , defined above, CoS0 o 


is given by 


Ir •> 


C05 &o ~ 


(23) 
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Since E q is independent of the incident energy, the energy 
dependence of Co S0 O is determined solely by ^ and . 

We have also assumed the same distribution function, 
given by equat ion(22) j for the multibody reactions npd,dd,pt 

3 

and r»He when pions are produced. 

3 

There are no measurements of d, t, or He production by 

o<* reactions. We have assumed arbitrarily that the breakup 

4 

products of each He nucleus effectively cluster together al- 
lowing us to treat the interaction as a two-body process with 
the two "clusters" having an isotropic angular distribution in 
the center of mass system. Thus, from equations (4) through (9)^, 
with m^ = m^ = m g = m r> we derive the distribution function for 
d. t. and He^ produced in<*e( reactions with and without pion 
production to be: 

},*,=,) = TTf (24) 

The measurements of Bailey (1956) of the spectra of 
deuterons and helium-3 produced by 190 Mev protons on carbon 
are shown in Figure (8^ ? normalized for unit integral. As can 
be seen these spectra have the same form as those shown in 
Figure ( 7 ) and can also be fitted by an exponential characteri- 
zed by 1/E = 0.14 (Mev/nucleon) -1 . Assuming that the triton 

distribution also has this form and that E Q is independent of 

3 

incident energy, the d, t, and He distributions for pCNO re- 
actions are the same as those given by equations (22) and (23). 
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P roduction Spectra 

The rate of production per gram of interstellar material 
of deuterons and helium-3 nuclei as a function of secondary 
kinetic energy per nucleon can now be calculated by using the 
demodulated cosmic ray spectra, production cross-sections and 
the secondary energy distributions discussed above. We have 
performed these calculations for a variety of modulating para- 
meters ^ and R q and in general the absolute production in- 
creases with increasing rj^ , but is quite insensitive to varia- 
tions in R q . The relative contribution of the various produc- 
tion reactions, however, is not strongly dependent on the as- 
sumed modulating parameters. For illustration, we show in 
Figures (9) and (lO) the deuteron and helium-3 production spectra, 
respectively, for = 350 Mv and R 0 = 0. As can be seen from 

Figure (9^, deuterons of energies less than about 70 Mev/nucleon 
are produced with equal probability by cosmic-ray helium-4 on 

4 o 

hydrogen in the pickup reaction H(He ,d)He and by cosmic ray 

protons in the multibody breakup of helium; deuterons in the 

energy range from about 70 to 200 Mev/nucleon are produced 
3 3 

mainly in the (-,-) resonance reaction of cosmic-ray protons 

2 2 

on hydrogen, H(p,d)tf* + ; and deuterons of energy greater than 
about 200 Mev/nucleon are produced with roughly equal probabi- 
lity by the multibody breakup reactions of cosmic-ray helium-4 
on hydrogen. It should also be noted that the«U( reactions, 
under the assumptions made above, contribute approximately 20% 
of the total deuteron production at all energies, and, because 
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of the lack of experimental data on the cross-sections and 
kinematics of these reactions, an uncertainty of this order 
is introduced in the calculation from this source alone. 

From Figure (lO) we also see that helium-3 nuclei of energies 
less than about 40 Mev/nucleon are produced principally by 
cosmic-ray helium-4 on hydrogen in the pickup reaction H(He 4 . 
He 3 )d, as are the low energy deuterons , while helium-3 nuclei 

of higher energies are produced predominantly by the breakup 

t 4 

of cosmic-ray helium-4 in the reaction H(He ,t)2p, the tritium 
decaying to helium-3. The eW reactions also produce about 20% 
of the helium-3 and introduce a corresponding uncertainty in 
the total production. 
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III. DEUTERON AND HELIUM-3 FLUXES AT THE EARTH 

From the deuterium and helium-3 production rates in the 
interstellar medium, calculated above, we may now determine the 
equilibrium density of these secondaries in the galaxy and then, 
taking into account modulation, the flux of the d and He 3 nuclei 
at the earth. The equilibrium densities of these nuclei in 
the galaxy can be obtained by solving a steady-state continuity 
equation in energy space which takes into account leakage from 
the galaxy, ionization losses and nuclear breakup of the secon- 
daries . 


h. + M 


Uj f u \j ■ ?V E) 


(25) 


where Ug is the secondary cosmic-ray equilibrium density: q is 

s 

the secondary production rate calculated above, per gram of 
interstellar material; ^ is the density of interstellar material 

O 

in g/cm ; ^ is the leakage lifetime from the galaxy; dE/dt 
is the rate of energy loss due to ionization; and A<j/v is the 
mean lifetime against nuclear breakup. The solution of this 
equation can be written as 


= lit! T»J«/*|] 


I je/Jtl 


(26) 


where T is the effective lifetime against leakage and nuclear 
breakup and is given by: 
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A more general expression for iig , which under an ap- 
propriate simplification reduces to equation (26) given above, 
can be obtained by introducing a probability distribution per 
unit time, Rt), such that PCt) dt is the probability that an 
observed particle (at t = G) was produced at a time t in the 
past. In terms of this probability, the number of particles 
of energy in dE around E which were produced in the time in- 
terval dt around t, with energies in dE* around E f , is given 
by 

(JU s O) = ^( E ‘) (28) 


where the energies E and E' , the time interval t and the 
energy loss rate dE/dt have to satisfy the following relation- 
ship: 




*5 

I 




I Jt/JM 


(29) 


E 


Equation (28) can be integrated by assuming that the rate of 
particle production $V E) is time independent. By changing 
the variable of integration from t to E' we find that the 
equilibrium density u g (E) is given by 
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f __ 




(30) 


Equation (30) is similar to equation (26) given above, ex- 
cept that the exponential distribution is replaced by the 
more general expression, P(t) . The functional form of this 
expression, however, depends on a variety of factors such as 
the spatial distribution of the sources, the nature of par- 
ticle propagation from the sources to the earth, and the 
properties of the trapping volume of the particles. 

In the present study we have assumed a simple exponential 
distribution, P(t) = e - */'* . Such a form is valid for the 
sudden losses resulting from nuclear breakup and it might be 
a good approximation to a physical situation in which the 
sources are uniformly distributed in a trapping volume of 
radius &^«mi which cosmic rays escape by diffusion caused by 
scattering off magnetic irregularities. The leakage lifetime 
, mentioned above, may then be written as ^ • «V(2^ 
where is the scattering mean free path. For such an ex- 

ponential distribution, equation (30) becomes equivalent to 
equation (26) . 

We now consider the numerical evaluation of equation (26). 
The energy loss rate dE/dt due to ionization can be replaced 
by BvdE/dR, where R(E) are charged particle ranges given by 
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Barkas and Berger (1964) . The deuteron breakup cross-section 
varies with energy but since we find that the lifetime Tj is de- 
termined principally by leakage, we have taken it to be a con- 
stant of 170 mb which corresponds to the average value over 
the energy range of 10 to 100 Mev/nucleon (Van Oers , 1963; 
and Davison et al., 1963). The helium-3 breakup cross section 
is much smaller yet and hence was neglected. The breakup mean 
path length, Aj > is given by , where N a is 

Avogadro’s number. 

Introducing a mean path lenght against leakage, X = J° 
equation (26) can be rewritten as 

co 

V 6) J At 'V Ww r[- 

c 




(31) 


where j (E) is the secondary cosmic-ray equilibrium intensity 
s 

in interstellar space. We have evaluated equation (31) for 

a variety of production spectra, q (e) , corresponding to a 

s 

range of values of the modulating parameters and R q , de- 

3 

fined above. We have then obtained the d and He intensities 
at the earth by remodulating the resultant equilibrium spectra 
in interstellar space with the same modulating function as that 
used in the demodulation process. The results of these calcu- 
lations, for R o =0 and for several values of the mean path length 
x and the modulating parameter Tj , are shown in Figures (11) 
and (12) for deuterons and helium-3 nuclei respectively. 


The net effect of the process of demodulation and the sub- 
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sequent remodulation, described above, is a decrease of the 
secondary intensity, because secondaries of a given energy 
are produced by primaries of higher energy which are there- 
fore less modulated than the secondaries which they produce. 

This can be seen in Figures (ll) and (l2) from a comparison of 
curves A and C, for zero modulation and modulation characterized 
by ^ = 0.35 BV. Also shown in Figures (11) and (12), are 
measurements at the earth of deuterons and He^ nuclei. In 
comparing these data with our calculations, we first consider 
the case of zero residual modulation. As can be seen from 
curve A, we find that for a cosmic-ray path length of about 
4g/cm and zero modulation, the calculated He intensities at 
energies greater than about 200 Mev/nucleon are in good agree- 
ment with the measured intensities (O’Dell et al, 1966; Biswas 
et al, 1967; Dennis et al , 1967). However, we also see that 

3 

at lower energies the calculated d and He intensities are much 
greater than the measured values (Fan et al, 1966a, b) , and we 
find that for zero modulation no constant value of x can give 
a fit to the data. This strongly suggests that a significant 
residual modulation exists at solar minimum. This conclusion 
is supported by the observed gradient in the cosmic-ray inten- 
sity between the orbits of the earth and Mars (O'Gallagher and 
Simpson, 1967; 0;Gallagher, 1967). 

The rigidity dependence of this residual modulation, however, 
is not well known and as discussed above, there are significant 
disagreements among the various experimenters as to the energy 
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dependence of the observed cosmic-ray spatial and temporal 
variations. According to Jokipii and Coleman (1968), the ob- 
served power spectrum of the interplanetary field variations 
indicate that the function f (R) , defined above, can be repre- 
sented by a power law in R with a varying spectral index, 
ranging from about -0.5 at low rigidities to -2 at high rigid- 
ities. Within the experimental errors, this would be consis- 
tent with the modulating function given by equation (2) with 
R q of the order of 600 MV. Such a rigidity dependent modu- 
lation is indeed demanded by the observed deuteron and helium-3 
data. As can be seen from curves A in Figures (11) and (12), 
the calculated d-to-He 3 ratio at, for example, 50 Mev/nucleon 
for zero residual modulation is about 1.8, whereas the measured 
ratio at the same energy per nucleon is perhaps as great as 6. 
Such a difference could result from a modulation of the form 
exp (”*^/Rp) > as given by equation (2) for R > R Q . Since at 
the same energy per nucleon, deuterons have a higher rigidity 
and hence are less modulated than He 3 nuclei, such a rigidity 
dependent modulation would increase the calculated d-to-He** 
ratio. Thus, the comparison of the d and He 3 measurements with 
curve A suggests not only the necessity of some residual modu- 
lation but also that above a rigidity of about 600 MV this mod- 
ulation is velocity as well as rigidity dependent. Because of 
this, we wish to suggest that the comparison of future measure- 
ments of the d-to-He 3 ratio with the calculations described in 
the present paper, could determine more accurately the charge 
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dependence and energy spectrum of the residual solar modulation. 

Since the present calculations do not depend on processes 
occurring below rigidities of about 500 MV we have calculated 
the d and He 3 spectra at the earth for R Q = 0 , and a range of 
values of x and i£, thereby determining the values of these 
parameters which best fit the observed intensities. Within 
the experimental uncertainties, we find that the available 
data on both the deuterons and helium-3 nuclei can be best 
fitted by an energy- independent path length, x = 4 + 1 g/cm , 
and a residual modulation of the form exp(- ?) characterized 
by 1^ = 0.35 + 0.15 BV. This can be seen by the best fitting 
curve C and the limiting curves B and D in Figures (11) and 
( 12 ) . 

The d/He and He^/He ratios derived from the present cal- 
culations are shown in Figures (13) and (14). Biswas et al 
(1966) have recently concluded that the energy dependence of 
the He^/He ratio can be accounted for only by an energy-depen- 
dent path length. However, as can be seen from this figure, 
our calculations show that both He^/He and d/He ratios can be 
explained by an energy- independent path length distribution 
coupled with a residual solar modulation, which above about 
600 MV is both velocity and rigidity dependent and is of the 
form exp(-^/Rji). 
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FIGURE CAPTIONS 


1. Deuteron production cross-section for the reaction 

H 1 (p,d)'TT + . 

2. Deuteron, triton and helium-3 production cross-sections 
for pjle 4 reactions. 

3. Deuteron, triton and helium-3 production cross-sections 
for pCNO reactions. 

4. Ranges of the Lorentz factors of the center of mass, and 
of the secondary particle in the center of mass^ y Q and 
respectively, which for the indicated reactions can con- 
tribute to given secondary laboratory Lorentz factors, ^ g . 
The solid lines define the ranges of and ^ s * which 
satisfy the inequalities (10) , and the dashed lines re- 
present ^ * as a function of 

5. Ranges of the incident kinetic energies per nucleon which 
can contribute to a given secondary energy for the indi- 
cated deuteron producing reactions. These ranges are de- 
termined from the requirements that cos 9 *£l and that 
the deuteron Lorentz factor in the center of mass be less 
than or equal to its maximum value determined from the con- 
servation of energy and momentum. 

6. Angular distribution in the center of mass of deuterons 
from the reaction He 4 (p,d)He 3 for various incident proton 
energies. 
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7. Energy distribution in the laboratory frame of deuterons 
and tritons from the multibody breakup of He 4 . 

8. Energy distributions in the laboratory frame of deuterons 
and helium-3 nuclei from proton-carbon interactions. 

9. Deuteron production spectra per gram of interstellar material 
for the various modes shown in Table I. These spectra were 
computed for a demodulated cosmic ray intensity with r|^=350 
MV and R o =0. 

10. Helium-3 production spectra per gram of interstellar material 
for the various modes shown in Table I. These spectra were 
computed for a demodulated cosmic ray intensity with ^=350 Ml/ 
and R q = O . 

11. The observed intensity of deuterons, together with the 
calculated deuteron intensities at the earth for various 
values of interstellar path lengths and modulating 
parameters . 

12. The observed intensity of helium-3 nuclei, together with 
the calculated helium-3 intensities at the earth for various 
values of interstellar path lengtlsand modulating parameters. 

13. The observed ratio of deuterons to helium nuclei, together 
with the calculated d-to-He ratios at the earth for various 
values of interstellar path lengths and modulating parameters. 

14. The observed ratio of helium-3 to helium nuclei, together with 
the calculated He-to^He ratios at the earth for various 
values of interstellar path lengths and modulating parameters. 



TABLE I 


DEUTERON PRODUCTION 

1 . p + H 1 -* d + ir + 

2. p + He 4 -*d + He 3 + (IT) 

3. o( + H 1 -»d + He 3 + (4T) 

4 . p + He 4 -* d + 2p+n + ( 4T ) 

5. <* + H 1 -*^ + 2p + n + (IT) 

6. o< + He 4 -*d + p + n +o( + (41') 

7. p + He 4 -*2d + p + (4T) 

8. o( + He^-*2d + p + (IT) 

9. o< + He 4 -* 2d +o(+ (IT) 

10. p + CN0-*d + ... + (If) 

11. CNO + H 1 -*^ + ... + (V) 

HELIUM-3 PRODUCTION 

1. p + He 4 -*He 3 + d + (IT) 

2. o< + H 1 -* He 3 + d + (IT) 

3. p + He 4 -*t + 2p + (4T) 

4 . ai + Hv — * t + 2p + ( if ) 

5 . o( + He 4 -* t + p+o( + ( IT ) 

6. p + He 4 -* He 3 + p + n + CtT ) 

7. o( + H 1 -* He 3 + p + n + (IT) 

8 . o( + He 4 -* He 3 + n + of + ( IT ) 

9. p + CNO -*He 3 +. . . + (H ) 

10 . CNO + H 1 He 3 . . . + ( IT ) 



TABLE II 


Kinetic Energy (Mev/nucleon) 


d producing reactions 

E th 

E z 

E min 

E s (th) 

E s (min) 

Eg(OP) 


pp-* d rr + 

286.5 

— 

312.4 

69.1 

63.9 

234.1 

p<* d He 3 

23.0 

27.7 

— 

0.9 

— 

230.2 

oi p -*d He 3 

23.0 

— 

37.3 

14.6 

9.1 

234.1 

p<X d 2pn 

31.7 

38.2 

— 

1.3 

— 

230.2 

o< p d 2pn 

31.7 

— 

51.6 

20.1 

12.5 

234.1 

po( -»dd p 

29.9 

36.1 

— 

1.2 

— 

230.2 

o( P -»dd p 

29.9 

— 

48.6 

19.0 

11.8 

234.1 

p O 16 -* d. . . 

14.3 

14.4 

— 

0.1 

— 

2847.9 

0 16 p-*d. . . 

14.3 

— 

27.1 

12.6 

6 . 6 

234.1 




He° producing reactions 

E th 

E z 

E min 

Eg (th) 

Eg (min) 

Eg(oo) 

po( -* He 3 d 

23.0 

37.3 

— 

0.9 

— 

37.8 

of p — v He 3 d 

23.0 

— 

27.7 

14.6 

12.2 

622.8 

pof-*t 

2p 

24.8 

40.2 

— 

1.0 

— 

37.8 

of p t 

2p 

24.8 

— 

29.9 

15.8 

13.1 

622.8 

pc* — > He 3 pn 

24.8 

40.2 

— 

1.0 

— 

37.8 

3 

oi P -*He pn 

24.8 

— 

29.9 

15.8 

13.1 

622.8 

p O 16 -* 

t . . . 

16.2 

16.4 

— 

0.1 

— 

1639.0 

1 

He 3 . . . 

» 

16.2 

16.4 

T“ 

0.1 

— 

1639.0 

o 16 P.l 

t . . . 

16.2 

— 

23.0 

14.3 

10.1 

622.8 

1 

'He 3 . . . 

16.2 

— 

23.0 

14.3 

10.1 

622.8 
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Figure 1 


KINETIC ENERGY (Bev) 



001 



(qui) -o 


Figure 2 


KINETIC ENERGY/ NUCLEON 




Figure 3 


KINETIC ENERGY/NUCLEON 











DEUTERON KINETIC ENERGY, E. (Mev/ Nucleon ) 



He 4 (n,n'i pt 

Tannenwald (1953) 
E = 90 Mev 

Innes (1957) 

E= 300 Mev 


exp {-0.27 E f } 


TRITON KINETIC ENERGY, E f (Mev/ Nucleon) 


Figure 7 






C(p,d) and C(p, He 5 ) 
Bailey (1956) E = I90 Mev 



Figure 8 


Deuteron and Helium-3 Kinetic Energy, E ( Mev/Nuc leon) 








DEUTERONS/M SEC. ST. MEV/NUCLEON 



DEUTERON KINETIC ENERGY (MEV/NUCLEON) 

Figure 11 



HELIUM-3 NUCLEI / M SEC.ST. MEV/NUCLEON 



HELIUM-3 KINETIC ENERGY ( MEV/NUCLEON ) 


Figure 12 



(9H/P) J 




r (He 3 /He) 





